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NOMENCLATURE 


C 

F(t) 

Fi(t) 

f,(t) 

K 

^1 

M 

«d 

«1 

u(t) 

Ud(t) 

Ui(t) 

Ur(t) 


Constant  matrix 
System  force  vector 
Disjointed  force  vector 
Element  force  vector 
Beam  force  vector 
Modal  force 

System  stiffness  matrix 
Disjointed  mass  matrix 
Element  stiffness  matrix 
Beam  stiffness  matrix 
System  mass  matrix 
Disjointed  mass  matrix 
Element  mass  matrix 
Beam  mass  matrix 
Nodal  displacement  vector 
Disjoint  displacement  vector 
Element  displacement  vector 
Beam  displacement  vector 
Modal  displacement 


Greek  Symbols 

6(t)  Delta  function 

6  Kroenecker  delta  function 

rs 

X  Eigenvalue 

X 

r 


Distinct  eigenvalue 


viii 


Natural  frequency 
Distinct  natural  frequency 
Eigenvector 
Distinct  eigenvector 
Normalized  eigenvector 
Inner  product 


1.  INTRODUCTION 


Recent  developments  In  the  control  of  flexible  structures  have 
created  the  need  to  experimentally  demonstrate  the  maneuver  of  flexi¬ 
ble  structures.  In  this  thesis,  an  experimental  set-up  for  a  flexible 
H-beam  type  structure  is  designed  (Figure  1). 

The  H-beam  model  is  designed  to  meet  several  objectives.  First, 
because  the  model  will  be  used  to  demonstrate  flexible  body  maneuver, 
the  H-beam  is  required  to  possess  significant  structural  flexibility. 
Secondly,  the  model  must  maneuver  about  three  degrees  of  freedom. 
Therefore,  the  structure  will  hang  from  a  cable  and  hinge  attachment 
at  the  center  of  gravity  of  the  H-beam.  Furthermore,  the  cable  must 
be  of  significant  length  in  order  to  make  the  pendulum  effect  on  the 
model  negligible.  Sizing  considerations  are  addressed  so  that  little, 
if  any,  static  deflection  (less  than  .01  inches)  or  torsional  deforma¬ 
tion  results. 
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2.  DEVELOPMENT  OF  EQUATIONS  OF  MOTION  OF  STRUCTURE 


The  equations  of  motion  for  the  H-beam  are  derived  using  the  finite 
element  method.  The  equations  of  motion  for  each  element  are  given  by 


«E 


Figure  2.  Finite  Element 


where  M^  and  K^,  are  the  element  mass  and  stiffness  matrices,  respec¬ 
tively,  for  axial  and  transverse  deformation  (Ref.  1).  For  axial  vi¬ 
bration,  linear  interpolation  functions  are  used.  Cubic  Interpolation 
functions  are  used  for  bending  vibration.  The  element  displacement  vec¬ 
tor  is  denoted  by  Ug(t),  and  the  element  forces  and  element  moments 
acting  at  the  2  nodes  of  the  element  are  denoted  by  F  (t).  The  element 
equations  of  motion  are  used  to  construct  the  equations  of  motion  for 
the  three  beams  comprising  the  H-beam  structure. 

The  equations  of  motion  for  each  beam  are  given  by 


M^  u^(t)  +  u^(t)  -  F^(t)  (1  -  1,2,3) 
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Figure  3.  Beams  Comprising  H-Structure  and  a 
Beam  Made  of  2  Finite  Elements 


where  and  are  the  mass  and  stiffness  matrices,  respectively,  for 
the  ith  beam,  and  F^(t)  and  u^(t)  denote  the  beam  force  vector  and  the 
beam  displacement  vector,  respectively.  The  disjointed  equations  of 
motion  for  the  structure.  Equations  (2),  can  be  express*  d  compactly 
as 

where : 


1. 


M.  0  0 

1 1  I 

- 1 - r - 

0  I  M,  I  0 

I  2  I 

of  o‘[«7 

I  I 


0  1  0 
.-_4- 


0  ^3 

I 


in  which  u^(t)  and  F^(t)  denote  the  disjoint  displacement  vector  and  the 


disjoint  force  vector,  respectively,  given  by 


Ui  (t) 


!l 


!d<^>  “  !2 


iJs 


h  <'> 


The  three  beam  displacement  vectors  comprise  Uj(t)  and  the  three  beam 


force  vectors  comprise  F,(t). 

«wU 


The  transition  from  the  local  coordinates  of  the  disjointed  system 


to  the  global  coordinates  of  the  structure  is  accomplished  using  the 


constraint  matrix  C  which  also  accounts  for  additional  constraints 


placed  on  the  structure  not  previously  considered.  In  particular,  it 


is  of  interest  to  satisfy  the  constraints  Imposed  at  the  boundaries  of 


the  elements  and  to  redefine  the  resulting  displacement  vector  in  global 


coordinates.  To  that  end,  the  disjoint  displacement  vector  is  related 


to  the  nodal  displacement  vector  by 


u  (t)  »  C  u(t) 

-•Q  -w 


where  u(t)  is  the  nodal  displacement  vector.  The  dimension  of  the  dis¬ 


placement  vector  is  identical  to  the  number  of  degrees  of  freedom  of 


the  system. 
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The  system  equations  of  motion  for  the  structure  are  obtained  from 
the  disjointed  equations  of  motion  by  substituting  Equation  (4)  into 
Equation  (3)  and  premultiplying  by  the  transpose  of  the  constraint  matrix 
to  obtain  the  system  equations  of  motion 

M  a(t)  +  K  u(t)  -  F(t)  (5) 

in  which 

M  =  M ,  C  ,  K  =  C  .  (6) 

d  d 

The  mass  matrix  M  and  the  stiffness  matrix  K  are  both  symmetric  n  x  n 
matrices  and  F(t)  denotes  the  system  force  vector 

F(t)  -  F.(t) 


In  assembling  the  equations  of  motion  for  the  H-beam,  it  is  assumed 
that  the  axial  deformation  along  the  three  beams  are  due  only  to  rigid- 
body  deflections  (Figure  4) . 

Associated  with  the  equations  of  motion  of  the  structure,  the 
eigenvalue  problem  is  defined  as 

XM<j)  -  .  (7) 


The  solution  to  the  eigenvalue  problem,  referred  to  as  the  eigensolution, 

is  comprised  of  the  eigenvector  (|)  and  its  associated  eigenvalue  X.  The 

matrices  used  in  this  problem  are  symmetric  n  x  n;  thus,  there  are  n 

linearly  Independent  eigenvectors  <j)^  (r  “  l,2,,,.,n)  and  n  associated 

eigenvalues  X^  (r  ■  l,2,...,n).  The  eigenvectors  are  often  called  the 

natural  modes  of  vibration  or  mode  shapes.  The  eigenvalues  are  related 

2 

to  the  structure's  natural  frequencies  .  by  X^  *  o)^  (r  »  1,2,.. .,n). 

I 
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The  lower  eigenvalues  can  be  computed  accurately  although  the  computed 
elgensolutions  with  higher  eigenvalues  are  grossly  In  error.  This 
error  Is  Inherent  In  the  discretization  process  that  takes  place  In  the 
modeling  of  the  structure  (Ref.  1).  Fortunately,  the  errors  are  not 
Important  provided  only  the  lower  modes  contribute  significantly  to  the 
system  response. 

The  numerical  solution  of  the  eigenvalue  program  Is  obtained  for 
the  eigenvalue  problem  In  the  form 


X<f  ”  A4> 


(8) 


Equation  (7)  must  be  manipulated  Into  the  form  of  Equation  (8).  Equa¬ 
tion  (7)  premultlpllng  by  M  the  Inverse  system  mass  matrix,  yields 
Equation  (8)  In  which  A  ■  M  ^K.  Associated  with  Equation  (8) ,  the 
orthonormallty  conditions  are  In  the  form 


(j)’’'’  <|)  -  6 

Xr  is  rs 


4)  A<|)  “  X  6 

Xr  -s  r  rs 


(9) 

(10) 


where  6  Is  the  Kronecker  delta  function  In  which  6  ■  0  for  t  ^  s 

rs  rs 

and  6^^  ■  1.  The  Interest  lies  In  normalizing  the  eigenvectors  with 
respect  to  the  mass  matrix  M,  such  that 


!r  ^ts 


rs 


(11) 


41'^  K4i  -  X  6 
Xr  -s  r  rs 


(12) 


s’, 


The  rigid-body  inodes  of  the  H  structure  are  linearly  Independent  and 
require  orthonormallzatlon.  To  that  end,  the  Gran>- Schmidt  orthogonallza- 
tion  process  Is  used  to  orthogonallze  the  eigenvectors  and  to  normalize 
the  eigenvectors  with  respect  to  the  mass  matrix  (Refs.  1,8). 

The  nodal  displacement  vector  u(t)  Is  a  linear  combination  of  the 
modes  and  written  In  the  form 

u(t)  -  I  <l>_  u_(t)  (13) 

r-l 

where  u^(t)  (r  ■  l,2,...,n)  are  modal  displacements  which  express  the 
extent  to  which  each  mode  contributes  to  the  overall  system  response. 

In  general  higher  modes  do  not  significantly  participate  in  the  re¬ 
sponse  so  they  can  be  excluded  from  Equation  (13).  The  modal  displace¬ 
ments  are  governed  by  the  scalar  equations 

u^(t)  +  J  u^(t)  -  f^(t)  (r  -  1,2 . s)  (14) 

where  s  <<  n  and  the  modal  forces  fj.(b)  are  related  to  the  nodal  force 
vector  F(t)  by 

f,(t)  -  4>3  F(t)  (r  -  1,2 . s)  .  (15) 

It  remains  to  compute  the  modal  displacements  in  Equation  (14). 

To  compute  these  modal  displacements,  the  rigid-body  responses  and 
flexible-body  responses  must  be  distinguished.  Because  impulse  forces 
are  the  type  of  modal  forces  used  on  the  structure,  new  initial  con¬ 
ditions  are  also  calculated. 
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2.1.  Rigid-Body  Modal  Responses  *  0) 


Equation  (14)  now  has  the  form 


u^(t)  -  *  1,2,. ..,ra) 


The  H-beam  structure  is  initially  at  rest  so  that 


u  (0)  -  0 
r 


u  (0)  =  0 
r 


The  nodal  force  vector  F(t)  is  expressed  in  the  form  of  an  impulse  as 


F(t)  -  F  6(t) 


The  rigid-body  response  is  computed  by  first  substituting  for  the  modal 
forces  in  Equation  (14) 


n  (t)  -  6(t) 

r  ~r  -o 


Introducing  Equation  (17)  into  Equations  (15)  and  (16)  and  Integrating 
over  the  time  of  the  impulse  from  0  to  0^,  yields 


rO  »  fO 

II  (t)  dt  -  (()  F  6(t)  dt 

Jo  Jo 


•  T 

u  (t)  *  (()  F 

r  g  ~r  ~o 


The  new  initial  conditions,  denoted  by  u  and  u  ,  become 

•'to  to' 


m  *4.  T 

u  u  (0  )  -  (p  F 
TO  r  Ir  ~o 


u  “  u  (0  )  *  0 
ro  r 


(18) 
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Equation  (16)  Is  now  replaced  with  the  homogeneous  system 


a^(t)  -  0 


Integrating  Equation  (19) ,  we  obtain 


(19) 


u^(t)  «  At  +  B 


u^(t)  =  A 


where  the  constants  A  and  B  are  found  using  the  new  Initial  conditions 
as  follows 


u  ■u(0)*0“B 
ro  r 


•  •  +  T 

u  =U(0)=(J)  F-A 
ro  r  ~r  ~o 


(20) 


The  rigid-body  response  Is 


(21) 


■  "ro‘  ■  '*2 . • 


2.2.  Flexible-Body  Modal  Responses  (o)^  4  0) 


The  new  Initial  conditions  due  to  Impulse  nodal  forces  are  computed 
using  the  same  Integration  scheme  that  found  u^^  and  u^^  for  the  rigid- 
body  responses.  Thus,  Equation  (14)  Is  Integrated  with  respect  to  time 


6(t)  dt 


0^ 

so  that 

+ 

rO 

1 

ii  (t)  dt  +0)^ 

r°  T 

u  (t)  dt  ■  4)  F 

Jo 

r  r  J 

0  r  -r  -o  J 

u  (t)  -  (()'  F 

r  -r  -o 


The  new  Initial  conditions  are  given  by 


•  •  +  T 

u  •  u  (0  )  *  (p  F 
ro  r  ~r  -o 


u  ■  u  (0  )  «  0 
ro  r 


Equation  (14)  is  replaced  with  the  homogeneous  system 


il  (t)  +0)  u  (t)  *  0 
r  r  r 


The  homogeneous  solution  has  the  form 


u  (t)  =  A  cos(u)  t)  +  B  sin((jJ  t) 
r  r  r  r  r 


u^(t)  =  -  A  0)  sln(a)  t)  +  B  uj  cos(a)  t) 
t  rr  r  rr  r 


where  constants  A  and  B  are  obtained  using  the  new  initial  conditions 
r  r  * 


u  *u(0)*0”A 
ro  r  r 


u  =  u  (O"*^)  •=  (})^  F  =  10  B 
ro  r  Ir  ~o  r  r 


B  ■  u  /u  -  ({)  F  /to 
r  ro  r  -r  ~o  r 


The  flexible-body  response  can  be  expressed  as 

u 

u  (t)  «  — ^  sin((jo  t)  (r  ■*  nri-1,  nri-2,  n) 

r 


(22) 


The  total  modal  response  of  the  H-beam  is  the  sum  of  the  rigid- 
body  responses  and  flexible-body  responses.  Introducing  Equations  (21) 
and  (22)  into  Equation  (13),  yields 


m  ^  n  u 

u(t)  ■y<j)ut+y(()  — ^  sin((o  t) 
~  -r  ro  _r,  -r  to  r 

r-1  nH-1  r 


(23) 


The  finite  element  code  created  models  the  H-beam  accurately  and 
the  accuracy  is  shown  by  comparing  the  natural  frequencies  of  a  six- 
element  model  with  a  twelve-element  model.  The  comparison  shows  the 
•  natural  frequencies  converging  and  the  Inclusion  principle  resulting. 

I 

)  To  test  the  nodal  response  accuracy,  the  six-element  model  is  excited 
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by  two  equal  and  opposite  impulse  forces  (.6  lbs  each)  applied  at  opposite 
ends  of  the  structure  (see  Figure  5).  The  response  at  the  various  nodes 
of  the  model  are  as  expected;  therefore,  knowing  that  there  is  a  good 
mathematical  model,  the  H-beam  can  now  be  dimenslonallzed. 


Figure  5.  H-Beam  Struck  with  2  Impulse  Forces 


Table  1.  Convergence  and  Inclusion  Principle  of  Natural  Frequencies 
for  6-Eleinent  and  12-Eleinent  Models 

6-Element  Model 

12-Element  Model 

6-Element  Model 

12-Element  Model 

6-Element  Model 

12-Element  Model 


612.36 

676.29 

1073.0 

1254.7 

650.2 

72A.8 

892.1 

1112 

1173 

1317 

1757 

82.113 

127.81 

301.85 

387.07 

570.78 

337.3 

359.1 

373.2 

396.5 

424.7 

469.2 

512.7 

554.6 

593.9 

0. 

0 

0.0 

0.0 

24 

.927 

43 

.278 

62.151 

0.0 

0.0 

0.0 

24.6 

42.8 

61.5 

80.1 

98.7 

121.1 

187.6 

292.6 

3.  SIZING  CONSIDERATIONS 


The  static  deflection  of  a  cantilever  beam  Is  considered  for  the 
sizing  of  the  beam  (Ref.  9).  For  a  specific  beam  composition,  aluminum 
for  example,  we  let  the  static  deflection  In  the  z  direction  Is  equal  to 
0.01  Inch,  and  the  weight  of  the  thrusters  Is  equal  to  eight  pounds  thus, 
the  beam's  length  can  now  be  calculated  for  a  specific  beam  cross  section 
(see  Table  2).  The  equation  below  is  used  to  compute  the  beam's  length 


Z 


ST 


(24) 


For  an  aluminum  beam,  the  following  values  are  given 


10  X  10^  psi 

^ST 

«  0.01  In 

.0975  Ib/ln^ 

P 

=  8  lbs 

pbh 

I 

12 

where  E  denotes  the  modulus  of  elasticity  (Young's  modulus)  for  aluminum, 
p  refers  to  the  density  of  aluminum,  Z^^  represents  the  static  deflection 
In  the  z  direction,  I  denotes  the  moment  of  Inertia  of  the  segment,  P 
represents  the  weight  of  the  thruster  assembly,  q  refers  to  the  weight 
per  unit  length  of  the  segment,  and  L  denotes  the  length  of  the  segment. 
The  length  and  cross  section  that  yield  the  best  combination  of  rigid- 
body  and  flexible-body  displacements  in  the  response  plot  are  the  design 
measurements  for  the  segment  (see  Figures  6,  7,  and  8).  The  beam  will 
have  the  same  cross  section  and  twice  the  length  of  the  segment.  Next, 
the  horizontal  beam,  connecting  the  two  parallel  beams,  must  be  sized. 


Table  2.  Parallel  Beams  Sizing  Table 
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A  hinge  assembly  that  connects  the  structure  to  the  hanging  cable 
from  the  celling  is  located  at  the  midpoint  of  the  horizontal  beam.  The 
cable  will  be  20  feet  long  to  negate  the  pendulum  motion  of  the  structure. 

j  Half  of  this  beam  is  taken  as  cantilevered.  These  two  segments  must  have 

I  maximum  rigidity  designed  in  the  z  direction  and  significant  structural 

flexibility  in  the  x  direction.  One  parallel  beam  and  two  thruster  assem- 

I 

j  biles  are  attached  to  the  free  end  of  these  segments.  The  static  deflec¬ 

tion  in  the  z  direction  is  set  at  0.01  inch.  For  ease  of  construction, 
these  segments  have  the  same  height  as  the  parallel  beams.  The  horizontal 

I 

j  segment’s  width  is  varied  to  find  the  length  and  cross  section  that  sat¬ 

isfy  the  static  deflection  criteria  (see  Table  3).  The  static  deflection 

;  is  computed  using  Equation  (24)  in  which 

P  -  2(8)  lbs  + 

I 

I  where  P  denotes  to  the  weight  of  the  two  thrusters  and  denotes  the 

weight  of  one  parallel  beam.  The  length  and  cross  section  yielding  the 
best  combination  of  rigid-body  and  flexible-body  modes  in  the  response 
plot  yield  the  design  dimensions  for  the  segment  (see  Figures  9  and  10) . 
The  horizontal  beam  is  twice  the  length  and  has  the  same  cross  section 
of  the  segment.  The  third  sizing  criteria  are  to  determine  if  torsional 
effects  enter  into  the  design  set  for  the  H-beam  structure. 

The  model  is  tested  to  see  if  torsional  modes  can  be  excited  by 
first  testing  to  find  the  static  angle  of  twist  on  the  parallel  beams. 

This  0g,j,  results  from  the  torque  created  by  the  weight  of  the  thruster 
assembly.  If  0^^  is  approximately  equal  to  zero,  the  static  angle  of 


Table  3.  Calculating  Length  of  Horizontal  Segment  for  2  Cross  Sections 


a.0 


twist  has  virtually  no  effect  on  the  structure.  The  static  deflection  for 
a  cantilever  beam  In  torsion  Is 


TL 

GJ 


(25) 


where  T  denotes  the  weight  of  a  thruster  assembly  times  the  maximum 
perpendicular  displacement  dj^  of  the  assembly  from  the  centroid  of  the 
beam's  cross  section,  L  represents  the  length  of  the  segment,  G  denotes 
the  beam's  shear  modulus,  and  J  denotes  the  moment  of  inertia  of  the 
cross  section.  Next,  the  fundamental  frequency  of  torsional  response 
Is  computed  to  rate  if  this  mode  might  be  excited. 

The  torsion  fundamental  frequency  is  calculated  for  the  desig¬ 
nated  parallel  beam  dimensions.  The  calculated  frequency  is  compared 
with  the  first  five  modal  frequencies  that  significantly  contribute  to 
the  nodal  response,  and  a  torsional  fundamental  frequency  significantly 
larger  than  these  modal  frequencies  indicates  that  the  structure  will  not 
be  Induced  Into  torsional  vibration  as  a  result  of  the  anticipated  struc¬ 
tural  excitations.  The  torsional  beam  fundamental  frequency  is  computed 
by 


where  3  equals  rTr/2L  for  a  cantilever  beam  and  I  denotes  the  mass  moment 
of  inertia  for  the  parallel  beam  cross  section.  The  H-beam  structure  is 
evaluated  to  see  if  this  torsional  mode  might  be  excited  (see  Table  4). 


Table  A.  Calculation  Torsional  0  and 
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4 .  HARDWARE  CON  S I DE  RAT I ON  S 


Next,  the  hardware  necessary  for  the  experimental  set-up  is  con¬ 
sidered.  Structural  dynamists  and  experimental  test  engineers  at  the 
NASA  Langley  Research  Center  in  Hampton,  Virginia,  are  consulted  and 
they  detail  several  requirements.  The  composition  of  the  H-beam  must 
be  specified.  In  addition,  certain  design  specifications  are  identified 
for  the  motors,  sensors  and  actuators,  and  hinge  assembly  used  as  experi¬ 
mental  equipment  for  this  set-up. 

4.1.  Composition  of  Structure 

An  H-beam  structure  composed  of  aluminum  is  recommended  because 
aluminum  is  light  weight,  one-third  the  weight  of  steel,  and  possesses 
excellent  material  properties.  Several  NASA  Langley  engineers  mentioned 
that  machinists  favor  aluminum  because  the  metal  is  "easy"  to  machine. 
Welding  properties  of  aluminum  alloy  6061  are  "first-rate."  Since  ease 
of  construction  is  paramount  in  building  the  H-beam,  aluminum  is  the 
choice  for  the  structural  composition. 

4.2.  Motors 

When  possible,  DC  permanent  magnet  motors  should  be  used  to  power 
the  sensors  and  actuators.  Because  these  DC  motors  have  only  two  wires, 
fewer  wiring  problems  are  Incurred.  Also,  DC  motors  are  easier  to  mount 
on  a  structure  than  AC  motors.  Motors  come  in  torque  drives,  brush 
drives,  and  permanent  magnet  drive  designs.  Permanent  magnet  DC  motors 


have  high  reliability,  long  operational  life,  and  require  little  mainten¬ 


ance.  To  power  the  air  thrusters  and  sensors,  24-volt  DC  motors  should 
be  used. 


4.3.  Actuators  or  Thrusters 

Thrusters  are  used  for  altitude  control.  Air  thrusters  are  utilized 
on  this  experimental  set-up.  A  thruster  assembly  is  mounted  at  each  end 
of  the  H-beam.  The  thrusters  are  positioned  to  fire  in  the  x-y  plane. 
Solenoid  actuated  cold  gas  jets  provide  the  impulse  forces.  A  compressed 
air  tank  supplies  air  to  the  jets  via  a  flexible  pressure  hose  running 
from  the  tank  up  to  the  ceiling,  down  the  steel  cable,  and  attached  to 
the  hinge  assembly.  The  air  passes  through  the  hinge  and  through  looped 
flexible  hoses  (Ref.  4)  down  to  each  thruster  assembly  which  consists  of 
four  solenoid  manifolds.  Thus,  the  air  leaves  the  flexible  hose,  goes 
through  a  four-way  valve,  and  through  short  hoses  to  each  of  the  mani¬ 
folds.  The  manifolds  gate  the  air  flow  between  the  regulated  supply  tank 
and  the  thrusters.  Thrust  is  initiated  with  a  discrete  command  that 
opens  the  solenoid.  These  thrusters  will  give  an  impulse  force  in  the 
range  of  .6  lbs.  On  the  next  three  pages  are  the  examples  of  the  type 
of  thruster  specifications  desired  (Tables  5,  6,  and  7)  (Ref.  2). 

4.4.  Sensors 

For  this  experimental  set-up,  three  types  of  sensors  are  required; 
position,  rate,  and  acceleration.  Position  sensors  detect  x  and  y  dis¬ 
placements  and  the  theta  angle  between  the  displacements.  Rate  sensors 


Table  5.  Thruster  Assembly  Specifications 


Thruster  System  (Example) 

Thruster  Manifold/Located  near  3-DOF  Hinge  with  3/8-inch  pipe  thread 
to  attach  air  hose 


Thruster  Assembly/Located  at  Nodes  1,  3,  4,  and  6 


Each  Assembly  Consists  Of : 


SPECS 


Thruster  Housing  -  4 
(see  schematics) 

Solenoid  Valve  and  Manifold  -  4 


4  -  1/4"  pipe  treading  (inlet) 

4  -  1/8"  pipe  treading  (outlet) 

Coil  rated  -  24  »/DC 


Nuraatics  Valve  Assembly  (on/off) 

(see  schematics  for  wiring  diagram 
for  solenoids) 


.25  Amps 
6  watts 

Operating  time  -  .028  sec 


Plumbing  Flared  Fittings  -  4 
Air  Hoses  -  5 

-  Parflex  PE 

-  Serial  #QC44J22C1 


Standard 
Wall  type  1 
Grade  4 
WP  125  psi 
3/8  OD  X  .062 


Valves  -  Rotary  two-position  type 

Usable  service  from  high  vacuum  to  300  psig  with  a  nominal 
cutoff  voltage  of  0.4V 
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Table  6.  Thruster  Data 


Typical  Thruster  Data  (Example) 


Supply  Pressure 

80  psi 

Peak  Thrust  (Impulse) 

0.6  lb 

Steady  State  Thrust 

0.3  lb 

Duration  of  Thrust 

24  sec 

Thruster  Control  Design  (Example) 
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detect  either  angular  or  displacement  velocities,  and  acceleration  sensors 


measure  linear  and  angular  acceleration  (Ref.  5). 


To  detect  linear  position,  rate,  and  acceleration,  servo- 


accelerometers  are  advised  for  this  set-up.  The  accelerometers  are 


mounted  as  close  to  the  thruster  assemblies  as  possible.  Rotational  rate 


sensing  units,  rate  gyros,  detect  angular  position,  rate,  and  accelera¬ 


tion  and  are  suggested  for  this  model.  All  sensors  used  on  the  model 


must  have  high  response  frequencies.  A  rate  gyro  attached  to  the  bottom 


of  the  H-beam  at  the  center  of  the  structure  will  sense  the  H-beam  going 


through  rotations.  On  the  next  two  pages  are  examples  of  recommended 


sensor  control  specifications  (Tables  8  and  9) (Ref.  2). 


4.5.  Hinge  Assembl3 


Attaching  the  structure  to  the  steel  cable  hanging  from  the  ceiling 


is  a  three  degree-of-f reedom  hinge  assembly.  This  assembly  must  easily 


support  the  weight  of  the  H-beam,  allow  the  structure  to  rotate  about 


all  three  axes  and  translate  in  any  two-dimensional  plane.  The  hinge 


must  also  allow  compressed  air  to  be  passed  through  it  to  the  four 


looping  flexing  hoses.  The  three  degree-of-freedom  movement  must  still 


be  possible.  Some  friction  may  result  in  the  hinge  because  of  the  air¬ 


tight  seals  required  in  the  assembly.  This  friction  will  Increase  the 


error  of  the  data  extracted  from  the  model.  A  great  deal  of  machining 


and  precision  engineering  is  necessary  for  such  a  complex  design. 


According  to  Jeffery  Williams,  a  NASA  structural  dynamlslst  who  fre¬ 


quently  uses  hinge  assemblies  in  his  control  of  flexible  structures 


experiments,  a  hinge  of  this  standard  costs  in  the  range  of  $5-10  thousand. 
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Table  9.  Rate  Gyro  Control  Specifications 
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5.  RESULTS 


After  comparing  the  numerous  test  cases  for  the  H-beam  dimensions, 
the  design  choice  is  an  H-beam  with  parallel  beams  four  inches  high,  3/8 
inch  thick,  and  70  inches  long,  and  a  horizontal  beam  four  Inches  high, 
3/8  inch  thick, and  58  inches  long.  This  design  shows  excellent  flexi¬ 
bility  and  easily  supports  the  anticipated  static  loads,  the  weight  of 
the  sensors  and  the  motors. 

Checking  the  possibility  of  torsional  excitation  on  the  design  pro¬ 
posed,  0  and  0)  are 
0 1  ^ 


..  (12  lbs)  (3/16  in)  (35  in)  „  „  _-6 

=  - 7 - X - 1—  =  9.8  x  10  rad 

^  (3.97  X  10°  Ib/in  )(2.02  in^) 


■>/(3. 97x10^  lb/in^)(12  in/ft)  (2.02  in^) 

1(1T) 

'i  1/12(. 004542  slug/in)(16.14  in^) 

[2(35  in)J 

5632.9  rad/s. 


The  static  angle  of  twist  is  approxiamately  equal  to  "zero"  for  this  design 
even  with  a  12-lb  thruster  assembly.  The  torsional  fundamental  frequency  is 
significantly  greater  than  even  the  fifth  natural  frequency  of  the  model. 
Therefore,  torsional  modes  will  not  be  induced  for  the  structure  proposed. 

Further  research  of  hinge  assemblies  concludes  that  the  criteria  set 
in  this  thesis  is  impractical.  As  previously  stated,  a  hinge  that  allows 
compressed  air  to  flow  through  it  requires  an  immense  quantity  of  engineer¬ 
ing  and  cost.  The  hinge  will  have  added  friction  resulting  from  the  air 
seals  on  the  hinge  and  will  hinder  motion  causing  a  significant  degree  of 
error  in  the  data  extracted  from  this  set-up.  For  cost  and  design  effec¬ 
tiveness,  a  simple  hinge  that  attaches  the  H-beam  to  the  steel  cable  and 
a  loosely  wound  flexible  air  hose  to  transport  compressed  air  from  the  air 


tank  to  the  hoses  on  the  structure  are  suggested.  NASA  aerodynamicists 
working  in  the  Langley  free-flight  wind  tunnel  state  the  stiffness  added 
to  the  model  because  of  the  compressed  air  in  the  hose  is  negligible. 
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6.  CONCLUSION 

A  finite  element  is  created  to  model  a  test  article  which  will 
demonstrate  control  of  flexible  structures.  The  H-beam  structure  will 
have  the  dimensions,  parallel  beams,  4"  x  3/8"  x  70";  horizontal  beam, 

4"  X  3/8"  X  70".  In  addition,  critical  hardware  specifications  are 
mentioned  for  the  motors,  thrusters,  and  sensors.  The  hinge  assembly 
proposed  allows  relatively  free  motion  with  errors  due  to  added  struc¬ 
tural  stiffness  being  negligible.  In  conclusion,  an  experimental  set¬ 
up  has  successfully  been  designed  to  demonstrate  flexible  body  maneuvers. 
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